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It is known that non-Abelian classical kinetic theory reproduces the Hard Thermal/Dense 
Loop (HTL/HDL) effective action of QCD, obtained after integrating out the hardest mo- 
mentum scales from the system, as well as the first higher dimensional operator beyond the 
HTL/HDL level. We discuss here its applicability at still higher orders, by comparing the 
exact classical effective action obtained in the static limit, with the 1-loop quantum effec- 
tive potential. We remark that while correct types of operators arise, the classical colour 
algebra reproduces correctly the prefactor of the 4-point function Tr only for matter in 
asymptotically high dimensional colour representations. 
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1. Introduction 



Most observables of QCD at a finite temperature T and chemical potential /i are not com- 
putable in perturbation theory beyond a certain order, due to severe infrared problems 
What can be done systematically however is the construction of various effective theories, 
obtained by integrating out only the hardest scales from the system. For static observables 
this leads to the concept of a dimensionally reduced effective theory Q, whose construction 
and non-perturbative properties have been studied in great detail (for reviews, see |3|). 

For non-static observables, on the other hand, the relevant effective description is the Hard 
Thermal and/or Dense Loop (HTL/HDL) theory |Q, ^ or, equivalently, its local reformula- 
tion 1^, 1^] as classical kinetic theory |9|. These constructions and their non-perturbative 
properties have however so far not been studied to a similar beyond-the-leading-order level 
as those of the dimensionally reduced theory (for a review of the current status, see [llO|| ). 

Very recently, there was a new positive indication on the effectiveness of the kinetic de- 
scription: it was noted that it reproduces also the first higher order operator beyond the 
HTL/HDL level [^], representing an effective (charge conjugation violating) three-gluon in- 
teraction. In this brief note, we wish to report on some properties of the classical kinetic 
theory at still higher orders. We point out that in this case agreement with quantum field 
theory is only obtained for matter in high dimensional colour representations. 

2. Formulation of classical kinetic theory 

We start by reviewing briefly the formulation of classical kinetic theory, used to describe how 
"hard" particles (quarks and gluons with momenta ~ max(T, //)), behave in the background 
of a "soft" (~ max{gT, gfi), where g is the gauge coupling) gauge field configuration A^. 

The starting point is to consider the hard modes as classical point particles carrying a 
colour charge Q", with dynamics governed by the Wong equations |^]. When the effect of 
collisions is neglected, the 1-particle distribution function obeys the Boltzmann equation 

Our sign conventions correspond to QCD with a covariant derivative in the fundamental 
representation = - igT^A"^, and in the adjoint Vf = d^'^d^ + gf^'^A''^; the field 
strength is Ff^^ = {i/ g)[D^,Dy\. We take Nf flavours of massless quarks and antiquarks, 
each carrying two helicities. For antiquarks one should replace Q —Q m Eq. (|^), and for 
gluons one should take the adjoint representation, Q — > Qadjoint = Q- 

For quarks/antiquarks the boundary condition for the solutions of Eq. (|^) is that for a 
vanishing gauge field background (say, at infinity), 

/(x,p,Q)^5+(p2)ny(u;pTM), <5+(p2) ^ 20(po) (2^) <5(p2) , (2) 
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where Uf is the Fermi-Dirac distribution function, and ujp = |p|. For gluons f{x,p,Q) — > 
6-^.{p^) ni){ujp), where is the Bose-Einstein distribution function. 

The solution of Eq. (|l|) defines a current induced by the coloured particles, 

jiix)=g J2 I PuQ'fix,P,Q), (3) 

helicitics "^P^Q 
species 

with 

dpQ f f f 



p J (27r) 7p 7p J (2-kY Jq J 
Here dQ is the colour measure, which contains delta functions fixing the representation de- 
pendent Casimirs (see, e.g., []T2|| ); we shall return to its properties presently. The current in 
Eq. @, in turn, defines via 

f.=-^fS„. (5) 

an action Sm = ^' + 6Sm, where = - /^(l/2)Tr F^^F^^ 

There are actually two somewhat different formulations of the kinetic theory, corresponding 
to the orders in which the integrals J^, Jq are to be carried out in the above. We shall mostly 
carry out first. However, if one takes the moments Jq {...), Jq Q"^ (...) of Eq. (||), imposing 
by hand the QCD type relation (say, for the fundamental representation) 

/ Q''Q'f{x,p,Q) = -Ls-bf{x,p) + Id'^'^rix^p), (6) 
Jq 2Nc 2 



with f = Jq f, = Jq Q""!, then one gets a closed set of equations for /, /" [| 

' dpv 



p-df-gp^F;^^^ = 0, (7) 



\ab fb 1 9 mbc^. 



We refer to this as the second formulation. Since the colour hierarchy is truncated by Eq. (^, 
the first and second formulations only agree for the low order colour moments. 



3. Exact solutions for special backgrounds 

Let us now recall that, for some special background field configurations, one can find an 
ansatz leading to an exact solution of Eq. (|^) . Such solutions have previously been discussed 
in the Abelian case [13, 14 1 and, for static backgrounds, in the non-Abelian H (see also ]l5[). 



However no comparison has been made with QCD, as far as we know. 

We will search for a solution of Eq. (||) in a form that only depends on the canonical 
momenta. Introducing the shorthands 

A'^g" = A, A^Q^ = A, (9) 
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we take for quarks (with a set of a's to be specified presently) 



f{x,p, Q) = + gA^). (10) 

Plugging into Eq. (|l]), we immediately find that if a's exist such that the condition daA^ = 
is satisfied for all a,v, then the ansatz in Eq. (|l^) solves Eq. (|l|). Similar solutions hold for 
antiquarks, by replacing Q ^ —Q, and for gluons, by replacing Q ^ Q. 

Apart from Eq. (||), the ansatz should also satisfy the boundary condition in Eq. (^. Thus, 
in the static limit {d^A'^ = 0), we obtain (for quarks) 

f{x,p, Q) = 6+{p^) nf{po + gAo - /u). (11) 

An exact solution could also be obtained in the homogeneous limit, diA'^ = 0, by replacing 
Po + gAo in Eq. (|ll|) by |p + ^A]. 

One may of course ask whether these exact solutions are the most general ones, given the 
boundary conditions in Eq. (|2|). It is at least easy to verify iteratively (without any ansatz) 
that they do agree with the perturbative solution in gQ^ of Eq. (|^) around f^^^ defined 
by Eq. (^): writing / = f^^^ + /^^^ + f^"^^ + the n^^ order solution is obtained by solving 

= , (12) 

where = 9^ + ^/'^'^Q'^^S^, L = gp'^Q'^F^.d^. 



4. Comparison with quantum field theory 

Let us now see what kind of an effective action, Sm, the solutions found in the previous 
section lead to, and compare with quantum field theory. 

Plugging Eq. (11) into Eq. (|3|) and taking into account the Nf flavours and the two helicities, 
we obtain, say, for the quark and antiquark contribution to the gauge current. 



Jo 



2gNf / nf{ujp + gAo - fi) - nf{u;p 

Jp,Q 

0. 



Solving now Eq. (^) and writing 5Sj^ = ^-C^^ , we obtain 



2NfT 



Inf 1 + (,^-'^P^9Ao+tJ-)/T\ ^{-wp+gAo-tJ.)/T 



The integral over p is easily carried out, and we finally arrive at 



(13) 
(14) 

(15) 



Hi 



Nf 



7 7r' 



1 



1 



2 90 127r2 



-.2 ,J'2 



g 43 



(16) 
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We have for completeness kept here even the field independent part, accounting for the leading 
1-loop ("free") expression for the fermionic contribution to the pressure of QCD. 

Going to Euclidean metric by writing Aq^ = iA^ , Ce = —^m{Aq^ — > iAQ), we observe 
immediately that Eq. (^) would agree with the fermionic contribution to the dimensionally 
reduced effective action [^] for Aq; or, equivalently, with the full 1-loop fermionic contribution 
to the effective potential for the phase of the Polyakov line ; provided that 



f Q'''Q''K..Q''" = hrT'^'T''^ ...T^"] , n = 0,...,4. (17) 

Jq L J symmetric part 

(The trivial factor (—1)" could be removed by inverting the sign convention for /i here, or 
in [0.) Thus the question is whether Eq. ([TtI) is satisfied. It is easy to show that it is for 
n < 3, while in general it is not for n = 4. 

Indeed, since the integration measure is gauge invariant |Q| and the Q^'s commute, the 
result of the left-hand-side of Eq. ( p!7| ) must have a covariant structure. For SU(2) this is of 
the form 

/ Q^'Q^^Q^^Q"^ = L(i?)(5"i"2^"3a4 ^ jaiaa ja2a4 ^ jaia4ja2a3-) _ 
JQ 

The integral here can be carried out even explicitly, for SU(2). Alternatively, to fix the 
constant L{R), we can contract this equation with 5"i°2^a3a4^ ^j^^^ g^j^ ^y^j- indices. Then 
the integral can easily be performed, due to the constraint 6{Q^Q'^ — C2{R)) where C2{R) 
is the quadratic Casimir, and the normalisation Jg = dn where d^ is the dimension of the 
representation. One finds L{R) = {dji/15)C2{R)- Therefore, in the classical effective action 
the piece quartic in Aq reads 

/ Al = \dRCl{R){A^,A^o)\ (19) 
while the quantum result, for an arbitrary representation of SU(2), can be seen to be 

TrA^, = ldRC2{R) {C2{R) - \] {AlAlf . (20) 



5 ^ ' V 3 

Writing da = 2j + 1, C2{R) = j {j + 1), we see that only for high values of j is the classical 
result in Eq. ([T9| ) a good approximation of the quantum one. (Suggestions along these lines 
were made already in ||8|, ^].) While the expressions given here hold for A'^c = 2, we expect 
similar conclusions for any SU(Ac), whereas for U(l) there is a perfect match. 

The reason why operators below the quartic one are correctly reproduced, is that the 
integral Jq{---) involves explicitly the quadratic and cubic Casimirs, which fix the symmetric 
parts of the traces of two or three SU(A^c) generators. The quartic term, on the other hand, 
depends on the commutation/anticommutation relations, and it is then at this order that 
one sees a difference between the quantum and classical colour algebras. 

While we have here focused our attention on the "first" formulation of classical kinetic 
theory, a similar conclusion can be reached for the "second" formulation, based on Eqs. (|7|), 
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(^), although the numerical discrepancy we find is different. In this case we do not know of 
an exact solution, but solve the equations iteratively, as was done to second order in |11]. At 
the third order, relevant for Tv Aq, we no longer find a solution local in space, unless we take 
^0 completely constant. Then, the result is too large by a factor three, for SU(2). 



5. Conclusions 

We have recalled in this brief note that for some special backgrounds, like a static one, the 
Boltzmann equation in Eq. (|l|) can be solved exactly, and the corresponding gauge field effec- 
tive action can be computed. We have then compared the result with the 1-loop dimensionally 
reduced effective action for quantum field theory. We remark that while the quadratic and 
cubic terms are correctly reproduced, the classical non-Abelian colour algebra generically fails 
at the next level. The relative discrepancy is the smaller the higher the dimensionality of the 
colour representation; however, for the physical QCD case, some more elaborate formulation 
seems to be required in order to have an exact match also at this level (see, e.g., [^]). 

Of course, all this does not imply that classical kinetic theory could not be a useful tool for 
a non-perturbative determination of many important observables of QCD, such as plasmon 
frequencies, damping rates, or physics related to the colour conductivity, for which effects from 
the quartic coupling are subdominant. Nevertheless, our observation should underline the 
need for a better beyond-the-leading-order understanding of the proper effective description 
of real-time dynamics in the QCD plasma, which would ideally at the same time also allow 
for a proper fixing of the renormalisation scale in the gauge coupling g as in the dimensionally 
reduced theory [|^, as well as for a systematic approach to the continuum limit [pO[ . 
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